Micellization of Sliding Polymer Surfactants 
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Following up a recent paper on grafted sliding polymer layers (Macromolecules 2005, 38, 1434- 
1441), we investigated the influence of the sliding degree of freedom on the self-assembly of sliding 
polymeric surfactants that can be obtained by complexation of polymers with cyclodextrins. In 
contrast to the micelles of quenched block copolymer surfactants, the free energy of micelles of 
sliding surfactants can have two minima: the first corresponding to small micelles with symmetric 
arm lengths, and the second corresponding to large micelles with asymmetric arm lengths. The 
relative sizes and concentrations of small and large micelles in the solution depend on the molecular 
parameters of the system. The appearance of small micelles drastically reduces the kinetic barrier 
signifying the fast formation of equilibrium micelles. 
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1. INTRODUCTION 

Predicting, controlling and finely tuning the self- 
assembly properties of amphiphiles through molecular 
design is a problem of central importance in physical 
chemistryi^ It is arguably also one of its major contribu- 
tions to other fields: in the biological realm, where self- 
assembled phospholipids build the walls of liposomes and 
cells; in cosmetics, pharmaceutics or detergency, where 
many formulations are self-assembled solutions of sur- 
factants, phospholipids and other amphiphile molecules. 
In this context, diblock copolymers have emerged as a 
paradigm for self-assembly. 3,4 ' 5,6 Typically, the insoluble 
block drives the chains to self-associate and the com- 
positional asymmetry between the soluble and insolu- 
ble blocks defines the assembling geometry* 5 - The many 
possibilities for architecture building offered by polymer 
synthesis and the development of polymer theory led to 
an unprecedented power of molecular control over the 
self- assembling structures of these so called macrosurfac- 
tants. It is nowadays possible, for instance, to build dif- 
ferent self-assembled structures from diblock, triblocks 
and many other block copolymers, to introduce charges 
at different places along the chains with single charge 
accuracy, to form reversible or frozen structures, to com- 
bine flexible and semi-flexible blocks, etc. However, a 
major difficulty still hinders progresses in the predictive 
power of micellization theories. Indeed, the macromolec- 
ular character of such surfactants implies that there is a 
large kinetic barrier for micelle formation, and the ther- 
modynamic predictions for typical quantities such as the 
critical micellar concentration or the aggregation num- 
ber are in many cases only marginally relevant. In this 
paper we discuss a novel macrosurfactant architecture, 
based on rotaxane inclusion complexes^ that can lead to 
a significant decrease of the kinetic barrier to micelliza- 
tion. 

Rotaxanes are molecular complexes formed when a 
ring like molecule, the rotor, is threaded over a linear 
molecule, the rotating axisi Unthreading of the ring can 
be prevented by subsequent capping of the axis ends. 8,9 



a) 





FIG. 1: Schematic representation of sliding a) unimers com- 
prised of soluble insoluble blocks of lengths N s and N c , re- 
spectively, b) small symmetric micelles and c) big asymmetric 
micelles. 



Rotaxanes can be made by combining different linear 
polymersi°*ii with different cyclic molecules, in several 
solvents^ One of the most well studied systems in- 
volve poly(cthylcnc-oxide) and a-cyclodextrins, which 
are oligosaccharides of 6 glucose units assembled as rings. 
When an inclusion complex formed by one a-cyclodextrin 
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and one PEO chain is attached to a surface by grafting 
the a-CD, it results a novel structure of polymer layers 
for which we recently^ coined the acronym SGP layers, 
for sliding grafted polymer layers. Contrary to the usual 
end-grafted polymer layer, in the SGP layers the chains 
retain the ability to slide through the grafting ring, a 
new degree of freedom that allows, as we recently have 
shown, to better relax the layer structure. In particular, 
our work suggests that in spherical, star-like layers, as 
obtained for instance if one attaches the a-CD ring to 
a small colloid, sliding of the polymers through the ring 
entails a truly versatile polymer layer structure. At low 
number of grafts the chains adopt mostly symmetric con- 
figurations with two comparable arm-sizes whereas with 
many grafts the adopted configurations are asymmetric 
with essentially one long arm participating in the corona. 
It was also shown that there exists an intermediate range 
of grafting numbers where the number of arms partici- 
pating in the corona is constant and where the addition of 
more grafts leads to the progressive replacement of sym- 
metric configurations by twice as much asymmetric ones. 
Conversely, in flat dense SGP layers, chains adopt only 
asymmetric configurations, having hence the same struc- 
ture as the quenched ones. An important consequence for 
self assembly is that the free energy of the curved SGP 
layers has a different structure from the usual polymer 
layers grafted to spherical surfaces, but that the free en- 
ergy of flat SGP layers only marginally differs from their 
non-sliding equivalents. 

We consider the micellization of decorated "one-pearl 
necklaces" where the driving force for self-assembly is 
provided by the bead decoration, an insoluble chain 
chemically attached to the ring molecule (see Figure^. 
The consequences of the sliding degree of freedom on 
the self-assembly of sliding polymeric surfactants are dis- 
cussed by focusing on curved self-assembled structures, 
since for flat self-assembled layers only marginal differ- 
ences are expected with the usual grafted layers. After 
a short reminder on diblock-copolymer micellization, we 
recall the theoretical description of the spherical SGP 
layer, before a detailed discussion of the possible micel- 
lization scenarios for this new family of macrosurfactant 
architectures (Figure ^) . 



2. MICELLIZATION OF 
DIBLOCK-COPOLYMERS 



In this section we briefly discuss the theory of micel- 
lization of block copolymers, see e.g. Ref. 14. Consider 
a solution of monodisperse block copolymer with N c in- 
soluble and N s soluble units. The insoluble blocks tend 
to aggregate and favor large aggregates while coronas of 
soluble blocks oppose the formation of big micelles. The 
free energy per unit volume F, of the solution of non- 
interacting micelles is a sum running over the aggregation 
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FIG. 2: Chemical potential fi as a function of total concen- 
tration for quenched (dashed line) and sliding (solid line) mi- 
celles. Parameters used: Ns = 300, N c = 62 and a = 0.5. 



number p: 



F = 
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where F p and c p is the free energy and the number density 
of a micelle comprising p surfactant chains respectively. 
The multiplier N T b 3 with N T = N c + N s being the total 
volume of a copolymer^ This multiplier was often over- 
looked, a general discussion is provided by Riessji& who 
gives the relevant elementary volume, in our case 7V<r& 3 , 
for microemulsions, assemblies of droplets, etc. The equi- 
librium distribution of chains in micelles is obtained by 
minimizing F with respect to c p along with the mass 
conservation constraint: 



(2) 



0=1 



where <p is the total number of chains per unit volume. 
This gives the equilibrium densities of p-arm micelles in 
the form c p b 3 NT — exp {(/xp — F p )/ksT}, where: 



H = F x + k B T\ii{c 1 N T b 3 ) 

is the chemical potential of unimers. 

Thus we can rewrite the equilibrium c p as 



(3) 



(c 1 N T b 3 ) p 
° p= N T b 3 6XP 

Rather than the distribution < 



(4) 



v > 



F p - pF\ 
k B T 

we will often use the 



more convenient function VL p = ln(ci/c p ): 

pFi - (p - 1) k B T In ( Cl N T b 3 ) (5) 



fi — F 
sup — ± P 

which may be considered as the thermodynamic poten- 
tial. Given a value of c\ in eq. 0] allows us to calcu- 
late the whole distribution c p and the corresponding to- 
tal chain concentration <fi. Experimentally, the chemical 



3 



potential /i is often studied as a function of <fi. We ob- 
tained the chemical potential curve from eqs.Q-Q with 
ci as the parameter. A kink in the representation (Fig- 
ure [3 dashed line) is commonly used as an indicator for 
the onset of micellization. Usually only a narrow range 
of unimer concentrations at equilibrium c\ falls in the 
physically acceptable condition: Nx^b 3 <C 1. When no 
interaction between micelles is accounted for, the more 
stringent criterion J2 p c p V p < 1 should be obeyed where 
V p is the volume of a micelle comprising p surfactants. 

The micellization scenario depends on the precise form 
of the free energy F p , discussed in the following para- 
graphs for quenched copolymers and annealed sliding 
copolymers. 



3. FREE ENERGY OF STAR-LIKE MICELLES 

The free energy of a micelle contains the molecular 
characteristics of the micellization process. Thus, in or- 
der to investigate the micellization of sliding polymer sur- 
factants, we have to specify the explicit form for the cor- 
responding F p . The free energy of a micelle is the sum 
of the core and of the corona contributions. Hereafter 
energies are expressed in ksT units. 



3.1. Core contribution 

We assume that the insoluble blocks form a dense ho- 
mogeneous core where soluble blocks and the solvent can- 
not penetrate. Hence the micelle free energy has two core 
contributions: (i) the surface tension term, F c = AiraR 2 , 
where R c is the radius of the core and a is the core-solvent 
surface tension expressed in ksT/b 2 units. For an incom- 
pressible core of size R c = (3/(4 7 r)piV c ) 1 / 3 6, this leads to 
F c = caN 2 ^ 3 p 2 ^ 3 , where c = (36-7T) 1 / 3 . (ii) Gaussian elas- 
tic contributions arising from stretching of the insoluble 
blocks in the core, F e i = wpR 2 c /{N c a 2 ) = wp 5 ^ 3 /N^ 3 , 
where w = 37r 2 /80 reflects the spherical geometry of the 
core. The value of w is obtained from SCF theory in the 
strong stretching limit. 17 Although the elastic contribu- 
tion of the core is only larger than UbT for micelles with 
large cores, we keep it for convenience^ 



expansions (e = 4 — d) are available, to first ordertiS 

lp-l = ~P(p-Z) + o(e 2 ) (6) 

Recently Monte Carlo simulations were carried out by 
Hsu et alm& in order to find the exact values of 7 P for a 
large range of p values. The results do not quantitatively 
agree with the classical predictions of the Daoud-Cotton 
modetfi 7 p — 1 ~ ~p 3 / 2 -)-..., indicating that the asymp- 
totic limit of very large p numbers, where subdominant 
powers of p are negligible, is not yet reached. If we insist 
on fitting Hsu et. al. results to the Daoud and Cotton 
model, say between p = 20 and p = 60, where the sim- 
ulated arms are still fairly long, the obtained amplitude 
is close to 0.2. It seems that the asymptotic limit, re- 
stricted to the leading term, is also out of experimental 
reach. Authors showed that the best fit of the simulation 
data with a power law ~ — p z is obtained with z = 1.68. 
Our discussion below is based on the simulation data and 
we will use, for convenience, the fitting function of the 
form close to eq. (0: —p(bp — a) z /16 reflecting the sim- 
ulation results. The best fit is 

7 P -l = -^(1.5p-6) - 7 , p>4 (7) 

while for p < 4 we take values from the table of Ref. |2jJ. 

When there is sliding degree of freedom, the free en- 
ergy of the corona has a more complex structure. The 
possibility for the soluble blocks to slide through a ring 
results in an annealed arm length distribution in the coro- 
nas of such micelles made of p-chains. Adjusting the 
length of the arm in the corona can decrease the crowd- 
ing effect originating from the steric repulsions between 
soluble blocks. In our previous papeui 3 - we found three 
possible regimes for the sliding coronas depending on the 
number of sliding chains per aggregate. The transition 
between these regimes is determined from a threshold 
value of the number of arms p* , defined by the conditions 
7p* — 7p*-i > — 1 and 7 P *+i — 7 P * < — 1. If the number 
of chains per aggregate is small, p < p*/2, the corona is 
fully annealed and the sliding chains arc likely to adopt 
any configuration. We call such micelles symmetric mi- 
celles, since symmetric configurations are prevalent. The 
corona properties are determined by configurations with 
2p arms and its partition function is given byi£ 



3.2. Corona contribution 

For large soluble blocks, as considered here, the corona 
is usually envisioned as a star of soluble blocks radially 
stretched away from the spherical core. The partition 
function Z p of a star with p equal arms of contour length 

N s is given by the critical exponent j p , Z p ~ N2 P 
Star exponents 7 P are known exactly in two dimensions 
and for ideal chains (d > 4). In three dimensions e - 



ZpOciVJiV?*- 1 , p<p*/2 (8) 

For a larger number of grafting chains, in the interme- 
diate regime, p* /2 < p < p* , only p* — p chains adopt 
symmetric configurations, while the rest of the chains are 
stretched out from the ring. The corona has p* arms and 
its partition function reads 

Z p oc Nf-PNy \ p*/2 <p<p* (9) 
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FIG. 3: Potential f2 p = ln(ci/c p ) as a function of the aggre- 
gation number p for sliding micelles (solid line, total surfac- 
tant concentration <f> — 3. x 10~ 8 ) and for quenched micelles 
(dashed line, <j> = 1. x 10" 10 , dash-dotted line, <f> — 3. x 10~ 8 , 
same as for sliding micelles) for the set of parameters: Ns = 
300, N c = 70 and a = 0.5. Micelles form almost without 
kinetic barrier from sliding polymer surfactants. 

For even larger micelles, p > p*, all the chains in the 
corona are strongly asymmetric and the partition func- 
tion of the corona coincide with the partition function of 
a quenched star: 

Z p <xN]*-\ p>p* (10) 

The determination of the threshold value of the num- 
ber of arms p* strongly depends on our knowledge of 
the values of the critical exponents j p . Our previous 
estimates^ were based on the first order e-expansion (Eq. 
|5J) yielding p* = 9. However, if we use the values obtained 
from the computer simulations in ref. l20l we get p* = 18. 
Thus, in the following we assume the latter estimate as 
most accurate. We use our fitting function Q for the 
values of j p , and the free energy of the sliding corona of 
a micelle Fp° rona = —\nZ p is determined as a piecewise 
function of the number of chains p. 

3.3. Total free energy 

Combining all terms together, the energy of a micelle 
is given by 

F p = F™ rona + ca7V c 2/ V /3 + wp^/N^ 3 (11) 

with 

F™" = f[p) In N s (12) 

where f(p) is defined piecewise by eas. (l8l9H0|) as dis- 
cussed in the previous section for sliding copolymers 
whereas for quenched copolymers: f(p) — — (7 P — 1). 

Note that the corona free energy of sliding surfactant 
unimers (p = 1) is F( = — 7ilniV s while that of the 



quenched unimer is = —(71 — 1) In _/V s . Thus the free 
energy of sliding surfactant is about — In N s less than 
that of quenched ones with the same soluble block size. 

Next we use the expressions of the free energy F p to 
compute the micelle size distribution for both quenched 
and sliding copolymers. 



4. MICELLE-SIZE-DISTRIBUTION: SLIDING 
VS. QUENCHED COPOLYMERS 

A preferred aggregation number is reflected by a min- 
imum in the potential £l p . Roughly speaking, the cor- 
responding aggregates will dominate over unimers if the 
minimum potential value is negative. For usual quenched 
copolymers the potential has only one, pretty sharp, min- 
imum (see Figure |3J and the polydispersity of spheri- 
cal micelles is usually very low, the distribution c p being 
sharply peaked around the average aggregation number 
p m - In this sense, quenched copolymer association should 
be termed closed association where the distribution is 
dominated by aggregates of uniform size with weak fluc- 
tuations appearing at a well defined unimer concentra- 
tion. In the common used classification, open association 
means that aggregates of different sizes significantly con- 
tribute to the distribution and appear gradually. The 
practical relevance of this classification has been criti- 
cally discussed recently^ 

The average aggregation number and the CMC are 
found from the conditions: f2 p = dVl p /dp = 0^ The 
minimum is separated from the origin by a kinetic bar- 
rier. Thus it takes a typical time ~ exp(— C/ m ax/^sT) 
for isolated chains to form a micelle; with J7 max the max- 
imum of the barrier. A rough scaling estimate gives 
C^max/fcT ~ aNc^ 5 ~ Depending on the conditions, this 
time for usual block copolymer micelles can be very large. 
Very often, quenched copolymer micelles do not form over 
reasonable time scales at the CMC, where they are ther- 
modynamically favored. Only at much higher concen- 
trations, the barrier U ma , x /kT becomes low (see Figure 

EJi. 

In contrast to the quenched case, the potential Cl p = 
ln(ci/cp) (Eq. JSJl) for sliding micelles can have two min- 
ima. One corresponds to small symmetric micelles with 
aggregation numbers up to 9, and the second corresponds 
to big asymmetric micelles (Figure yj. The appearance 
of small micelles leads to a drastic decrease of the kinetic 
barrier, such that the first minimum is separated from 
the origin by a barrier of order kT signifying the fast 
formation of equilibrium micelles. The chemical poten- 
tial \x of sliding copolymers as a function of (f> presents a 
rounded kink as compared to the quenched case (Figure 

Depending on molecular parameters of the system, the 
formation of small micelles can follow (Figure^ and^J;) 
or precede (Figure ^Js and |1JI) the formation of big mi- 
celles. Typical examples of the mass distribution in mi- 
celles, pc p /4>, with increasing polymer concentration are 
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FIG. 4: Typical examples of the potential fl p = ln(ci/c p ) 
for different set of parameters: a) Ns = 300, N c = 60, a = 
0.5, <t> = 3.3 x 10~ 5 , b) N s = 300, iV c = 55, a = 0.6, $ = 
2.0 x 10" 8 , c) N s = 1000, iV c = 60, a = 0.5, (j> = 2.8 x 
10" 5 , d) iV s = 10000, iV c = 60, a = 0.7, 4> = 7.8 x 10" 7 . 
Depending on parameters, small micelles appear first, a) and 
c), or formation of big micelles happens earlier, b) and d). 

shown in Figure [SJ The equilibrium aggregation num- 
ber switches between symmetric (small) and asymmetric 
(big) micelles as the concentration increases. Tuning the 
parameters of the system, either the surface tension of 
the core, a, or the relative length of the blocks (in our 
case we vary N c keeping N s constant), we can change 
the order of the appearance of small and big micelles. 
Figure IH1 shows the corresponding state diagram which is 
related to Figure [7| When big and small micelles coexist 
the big ones usually dominate by mass. Hence we rep- 
resent the boundary between big and big+small micelles 
by a dashed line. 

Though the full potential (Eq|SJl is uneasy to handle 
analytically, we may somewhat simplify it by omitting 
the stretching of the insoluble blocks. The aggregation 
number of a vanishing minimum of this approximate po- 
tential (for some chemical potential) is then solution of 
the following equation: 

(p-l) dJ ^-f{p) + .f(l) _l caNf 3 
p-z/3 + 2p -i/3 _ 3 3 In 7V S 

the left hand side of this equation h(p) collects all the 
terms depending on p. Thus, in this approximation the 
desired aggregation numbers depend only on the combi- 
nation of parameters (3 = ^ccrN^ 3 / In N s . It is easy to 
see that in the quenched copolymer case there is indeed 
only one solution of p m and in the Daoud-Cotton limit 

/ 2 / 3 \6/5 

Pm ~ ( o-N c /mN s J , a classical result. 

It is worthwhile discussing the number of solutions of 
cq. Ijl 3(1 as a function of the parameter (3 for sliding 
copolymers. It can be obtained by counting the number 
of intersections of h(p) (Figure 0) with the line parallel 
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FIG. 5: Evolution of the mass distribution in micelles, pC p /(f>, 
with increasing total concentration cj>- a ) 2. x 10 -10 , b) 5. x 
10" 10 , c) 2. x 10 -9 , d) 8. x 10 -8 for 7V S = 300, iV c = 70 and 
a = 0.5. 



to the abscissa. In the first region, (3 < 1.38, there is 
only one intersection. It corresponds to one minimum 
in £l p (p). The number of arms p is small, thus this re- 
gion corresponds to the coexistence of unimers with small 
symmetric micelles. In the region 1.38 < (3 < 2.77 there 
are three intersections showing the existence of two min- 
ima separated by a barrier. Small micelles with aggrega- 
tion numbers p < p* /2 = 9 coexist with big micelles with 
p > p* — 18. In the third regime, (3 > 2.77, there are 
only asymmetric micelles of high aggregation numbers. 
The crossover values of the parameter (3 are compati- 
ble with Figure where the crossover values of 7V C are 
40 and 93. For given values of parameters (N c , N s , a) 
it should further be checked that the micelles exist for 
physical concentrations. 

Figure [7| also shows the aggregation number of 
quenched copolymer micelles as a function of (3 (dashed 
line). The difference with large sliding micelles is due 
to a shift in the free energy of a unimer from F( to F*. 
The Daoud-Cotton line is also shown (dotted line). As 
already discussed, no agreement can be obtained using a 
simple Daoud-Cotton power law. 
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FIG. 6: Schematic diagram of different types of micelles 
present in a solution. Labels designate regions where the indi- 
cated structures are dominant. Parameters used: Ns = 300, 
a = 0.5. The boundary between the region of big micelles and 
the region of the coexistence of small and big micelles is in- 
dicated as dashed line, because big micelles usually dominate 
by mass in both regions. 

5. CONCLUDING REMARKS 

Sliding-diblock-copolymer surfactants where the solu- 
ble and insoluble blocks are topologically tethered by a 
small ring polymer show a much readier micellization be- 
havior than the corresponding quenched diblock copoly- 
mers. When the size of insoluble blocks is not too long, 
both small micelles (with an aggregation number slightly 
smaller thanp*/2 = 9) and large micelles (with an aggre- 
gation number slightly larger than p* = 18) may coexist 
in the solution. In the small micelle each copolymer par- 
ticipates in the corona with two arms whereas in the large 
micelle the asymmetric arm length distribution is dom- 
inant. Hence coronas in small and large micelles com- 
prise a similar number of arms. Small micelles can form 
without appreciable kinetic barrier at the lowest concen- 
tration where they are thermodynamically favored. The 
barrier opposing the formation of the larger micelle re- 
mains modest. This is in marked contrast to the cor- 
responding quenched copolymers where kinetic barriers 
are very high at the same concentrations. Which micelle 
appears first with increasing concentration depends on 
the block asymmetry. At somewhat higher concentra- 
tions, the micelle-size-distribution broadens and covers 
the whole range from the small to the large aggregation 
number described previously. 

For large insoluble blocks only the large micelles, sim- 
ilar to the quenched ones, form. Due to the loss of the 
sliding degree of freedom upon association, the annealed 



CMC is higher than the quenched one. On the other 
hand the kinetic barrier is markedly lower for annealed 
copolymer micelles (typically by a factor 2) that may 
form at the CMC within experimental times. 
4.0 n 

h (P) 3 .5: 




p 

FIG. 7: Graphical solution of the eq. (1 1^1 . For values of 
/3 = |ccrAc' /3 / In N s between 1.38 and 2.77 a small symmetric 
micelle and a large asymmetric one can be more stable with 
respect to the unimer (for different chemical potentials). The 
bold dashed line (p < 3) was not calculated. The thin dashed 
line corresponds to quenched copolymers, where there is only 
one micelle size. The dotted line corresponds to the Daoud- 
Cotton model for quenched copolymers. 



Our description of star-like micelles is based on excel- 
lent values of the vertex exponents from a recent simu- 
lation by Hsu et. alM We think that this improves over 
the standard Daoud and Cotton model for the classic 
quenched micelles. 

Sliding copolymers are a new interesting example 
of non ionic annealed copolymers. Some charged 
micelles can also display a similar behavior as re- 
ported by Zhulina and Borisov for insoluble/ annealed- 
polyelectrolyte diblocks^ To our knowledge only flat 
brushes of complexing polymers have been studied 
theoretically. 26 Diblocks where the soluble block forms 
a complex with small colloids (proteins) also belong to 
the class of annealed copolymers. One may speculate 
whether such diblocks also form two types of micelles. 
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